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We present estimates of light quarks masses using lattice data. Our main results are based on a global analysis 
of all the published data for Wilson and Staggered fermions, both in the quenched approximation and with 
rif — 2 dynamical flavors. The Wilson and Staggered results agree after extrapolation to the continuum limit for 
both the n/ =0,2 theories. Our best estimates, in the MS scheme at scale 2 GeV, are m = 3.2(4) MeV and 
m s — 90(20) MeV in the quenched approximation, and m ~ 2.7 MeV and m s ~ 70 MeV for the 71/ = 2 theory. 
These estimates are significantly smaller than phenomenological estimates based on sum rules, but maintain the 
ratios predicted by chiral perturbation theory (%PT). Along with the new estimates of 4-fermion operators, lower 
quark masses have a significant impact on the extraction of e'/e from the Standard Model. 



1. LIGHT QUARK MASSES 

The masses of light quarks m u , rrid, and m s 
are three of the least well known parameters of 
the Standard Model. These quark masses have to 
be inferred from the masses of low lying hadrons. 
%PT relates the masses of pseudoscalar mesons 
to m u , md, and m s , however, the presence of the 
unknown scale \i in £-^pT implies that only ratios 
of quark masses can be determined. For example 
2m s /(m u + TTid) = m s /m — 25 at lowest order, 
and 31 at next order [HQ. Latest estimates us- 
ing QCD sum rules give m u + md — 12(1) MeV 
|||. However, as discussed in Q], a reanalysis of 
sum rules shows that far more experimental in- 
formation on the hadronic spectral function is 
needed before sum rules can give reliable esti- 
mates. Thus, lattice QCD is currently the most 
promising approach. 

To extract a, m, m s we fit the global data as 



M PS = B PS {m 1 +m 2 )/2 
My = Ay+B v (m 1 +m 2 )/2 



(1) 



for each value of the lattice parameters, f3, rif, 
fcrmion action. From Bps, Ay, By we determine 
the three desired quantities; the lattice scale a 
using M p , m using M 2 /M 2 , and m s in three dif- 
ferent ways using Mk, Mk*, M$. Throughout 
the analysis we assume that <f> is a pure ss state. 
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Note that using Eq. |l| means that we can predict 
only one independent quark mass from the pseu- 
doscalar data, which we choose to be m. The 
reason for this truncation is that in most cases 
the data for M„ and M p exist at only 2 — 4 values 
of "light" quark masses in the range 0.3m s — 2m s . 
In this restricted range of quark masses the ex- 
isting data do not show any significant deviation 
from linearity. One thus has to use My in order 
to extract m s . Details of our analysis and of the 
global data used are given in ||. 

For Wilson fermions the lattice quark mass, 
defined at scale q* , is taken to be rriL(q*)a = 
(1/2k— l/2/c c ). For staggered fermions m^q*) = 
mo, the input mass. The MS mass at scale /i is 
mjj-^(p) — Z m (jia)mi,(a), where Z m is the mass 
renormalization constant relating the lattice and 
the continuum regularization schemes at scale \x, 
and A = g 2 /16tt 2 '. In calculating Z m , a la Lepage- 
Mackenzie, we use otjjs for the lattice coupling, 
use "horizontal" matching, i.e. \i = q* = 1/a, 
and do tadpole subtraction. We find that the 
results are insensitive to the choice of q* in the 
range 0.86/a — ir/a and to whether or not tad- 
pole subtraction is done. Once rn-jjg^) has been 
calculated, its value at any other scale Q is given 
by the two loop running. We quote all results at 
Q = 2 GeV. 

We extrapolate the lattice masses to a = us- 
ing the lowest order corrections (Wilson are 0(a) 
and Staggered are 0(a 2 )). In the quenched fits 
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Table 1 

Summary of results in MeV in MS scheme at 
H = 2 GeV. The label W(0) stands for Wilson 
with rif = 0. An additional uncertainty of ~ 10% 
due to the uncertainty in the lattice scale a is 
suppressed. 





m 


m 8 (M K ) 


m s (M ) 


m s {M* K ) 


W(0) 
S(0) 


3.3(4) 
3.1(1) 


83(10) 
78 ( 3) 


96(10) 
96 ( 2) 


76(20) 
87( 2) 


W{2) 
S{2) 


2.5(3) 
2.9(3) 


63(8) 
73(8) 


77(10) 
66( 6) 


78(22) 
59( 6) 



we omit points at the stronger couplings (a > 
0.5 GeV -1 ) because we use only the leading cor- 
rection in the extrapolation to a = 0, and because 
the perturbative matching becomes less reliable 
as is decreased. The bottom line is that we 
find that the leading corrections give a good fit 
to the data, and in the a = limit the two differ- 
ent fermion formulations give consistent results. 
Our final results are summarized in Table [j]. 

The global data for m and the extrapolations 
to a — for Wilson are shown in Fig. [|. Using the 
average of quenched estimates given in Table we 
get 

m(MS, 2 GeV) = 3.2(4) (3) MeV (quenched). (2) 

where the first error estimate is the larger of the 
two extrapolation errors, and the second is that 
due to the uncertainty in the scale a. 

The pattern of O(a) corrections in the un- 
quenched data (n/ = 2) is not clear and we only 
consider data for > 5.4. The strongest state- 
ment we can make is qualitative; at any given 
value of the lattice spacing, the n/ = 2 data lies 
below the quenched result. Taking the existing 
data at face value, we find that the average of 
the Wilson and staggered values are the same for 
the choices > 5.4, > 5.5, or > 5.6. We 
therefore take this average 

m(2 GeV) « 2.7 MeV (n f = 2 flavors) , (3) 

as the current estimate. To obtain a value in the 
physical case of rif = 3, the best we can do is 
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Figure 1. m(MS, 2 GeV) extracted using M, 
data with the scale set by M p . 



to assume a behavior linear in n/ . In which case 
extrapolating the n/ = and 2 data gives 

m(2 GeV) « 2.5 MeV (n f = 3 flavors). (4) 

We stress that this extrapolation in nf is ex- 
tremely preliminary. 

We determine m s using the three different 
mass-ratios, M^/M^, M K */M p , and M^/M p . 
Using a linear fit to the pseudo-scalar data con- 
strains m s (MK) = 25m. Using the vector mesons 
Mfr and M$ gives independent estimates. The 
quenched data and extrapolation to a — of 
m s (M < p) are shown in Fig. The average 

of Wilson and staggered values are m s (M < f > ) = 
96(10) MeV and m s {M K >) = 82(20) MeV where 
the errors are taken to be the larger of Wil- 
son/staggered data. From these we get our final 
estimate 

rn s = 90(15) MeV (quenched) . (5) 

The rif = 2 data shows a pattern similar to that 
for to. Therefore, we again take the average of 
values quoted in Table |l] to get 

to s = 70(11) MeV (n/=2). (6) 
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Figure 2. Comparison of m s (MS, 2 GeV) ex- 
tracted using Mff, for the quenched Wilson and 
staggered theories. 



The error estimate reflects the spread in the data. 

Qualitatively, the data show three consistent 
patterns. First, agreement between Wilson and 
Staggered values. Second, for a given value of 
a the rif = 2 results are smaller than those in 
the quenched approximation. Lastly, the ratio 
m/m s (M^) is in good agreement with the next- 
to-leading-order predictions of chiral perturba- 
tion theory for both the n/ = and 2 estimates. 
It is obvious that more lattice data are needed to 
resolve the behavior of the unquenched results. 
However, the surprise of this analysis is that both 
the quenched and rif = 2 values are small and lie 
at the very bottom of the range predicted by phe- 
nomenological analyses 

2. CP VIOLATION and e'/e 

A detailed analysis of 4-fermion matrix ele- 
ments with quenched Wilson fermions at = 6.0 
is presented in |6). The methodology, based on 
the expansion of the matrix elements in powers 
of the quark mass and momentum, is discussed 
in Q. Our estimates in the NDR scheme at 
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0.68(4) , 

0.78(1) , 

0.58(2) , 

0.81(3) . 



(7) 



The errors quoted are statistical. The major re- 
maining sources of errors in these estimates are 
lattice discretization and quenching. 

To exhibit the dependence of the Standard 
Model (SM) prediction of e'/e on the light quark 
masses and the B parameters we write 



e'/e = Al c + c 6 B 6 1/2 M r + c 8 B^M, 



? 3 / 2 / 



(8) 



where M r = (158 MeV/(m s + rrid)) 2 . For rea- 
sonable choices of SM parameters Buras et al. 
estimate A = 1.3 x 10~ 4 , Co = —1.3, cq = 7.9, 
C8 = —4.3 B. Thus, to a good approximation 
e'/e oc M r ; and increases as B^ 2 decreases. As a 
result, our estimates of m, m s ,B^ 2 increase e'/e 
by roughly a factor of three compared to previous 
analysis, i.e. from 3.6 x 10~ 4 to ~ 10.4 x 10~ 4 . 
This revised estimate lies in between the Fer- 
milab E731 (7.4(5.9) x 10" 4 ) and CERN NA31 
(23(7) x 10~ 4 ) measurements and provides a sce- 
nario in which direct CP violation can be ex- 
plained within the Standard Model. 
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